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Quantum phase transition of bosons in a shaken optical lattice 
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Recently, the lattice shaking technique has been used to couple different Bloch bands resonantly. 
For the one-dimensional (ID) case, in which shaking is along only one direction, experimental 
observation of domain-wall formation has been explained by superfluid Ising transition. Inspired by 
these, we generalize to a 2D case in which shaking is along two orthogonal directions. Analogous to 
the ID case, we find three different phases, the normal superfluid (NSF) phase, the D 4 symmetry¬ 
breaking superfluid (D4SF) phase and the Mott insulator (MI) phase. Furthermore, we demonstrate 
that the interaction effect induced by inhomogeneous band mixing can modify the critical shaking 
amplitude. Unlike in the ID case, shaking types also can modify the critical shaking amplitude. 
Unlike in the ID case, shaking types also can modify the critical shaking amplitude. We also 
construct a low-energy effective field theory to study the quantum criticality of bosons near the 
tricritical point of NSF, D4SF and MI phases. Moreover, we find a Bose liquid with anisotropically 
algebraic order and propose to change the Bose-Einstein condensation (BEC) into a non-condensed 
Bose liquid by tuning the shaking amplitude approaching the critical value. 


I. INTRODUCTION 

More and more interest has been attracted to ultra¬ 
cold atoms trapped in a time-periodically driven opti¬ 
cal lattice. There are two cases, off-resonance and res¬ 
onance. For the off-resonance case, it was stated that 
shaken lattice system can be described by an effective 
time-independent Hamiltonian with renormalized hop¬ 
ping amplitudes for a large shaking frequency [I]. It was 
also experimentally demonstrated that hopping ampli¬ 
tude can be changed dynamically with maintained phase 
coherence of condensation by shaking the lattices [2]. The 
lattice shaking technique can be used to tune hopping 
parameters and even invert the signs in a coherent way, 
which opens a new direction to simulate quantum phase 
transitions in ultracold atom systems. Coherent control 
of the superfluid-Mott-insulator (-MI) phase transition 
has been realized in a shaken three-dimensional optical 
lattice [3]. A synthetic gauge field can be realized in a 
shaken optical lattice 0H3, and this is equivalent to inset¬ 
ting a 7 r flux in each plaquette in a shaken square lattice, 
which generates a staggered-vortex superluid state [4], 
or in each triangle in a shaken triangular lattice, which 
generates various types of frustrated states 0E|. In ad¬ 
dition, interparticle interaction can be tuned from repul¬ 
sive to attractive in fermionic lattice systems by ac forc¬ 
ing, which allows one to simulate an attractive Hubbard 
model effectively with temperatures below the supercon¬ 
ducting transition temperature [8]. 

The resonance case starts from experimental obser¬ 
vation of domain-wall formation for bosons condensed 
in a shaken one-dimensional (ID) optical lattice [9], in 
which the lattice shaking technique hybrids different 
Bloch bands. The effective Hamiltonian cannot be de¬ 
scribed by renormalized hopping amplitudes or interac¬ 
tions as used in off-resonant cases, which may lead to 
novel phases. The finite-momentum superfluid phase 
with spontaneously broken Z 2 symmetry called the Z 2 SF 
phase has been observed [9] and the corresponding nor¬ 
mal superfluid- (NSF-)Z 2 SF-MI phase transition has 


been described by a low-energy effective field theory in 
the ID case m- 

The finite-momentum condensate has been realized by 
spin-orbit (SO) coupling generated by Raman transitions 
mm , or in a staggered magnetic field [T3], or in a 
shaken optical lattice BUM- The condensate with fi¬ 
nite momentum has spatially inhomogeneous order pa¬ 
rameter, which is a bosonic analogy to the Fulde-Ferrell- 
Larkin-Ovchinnikov phase in superconductors [13] . In¬ 
spired by the discovery of a finite-momentum condensate 
by resonantly shaking a lattice along one direction El HD], 
in this paper, we generalize to a two-dimensional (2D) 
case. Using Floquet theory, we demonstrate formally and 
numerically that lattice shaking leads to a phase tran¬ 
sition from the NSF phase to the D 4 SF phase as the 
shaking amplitude increases. We further show inhomo¬ 
geneous band-mixing-induced interaction effect modifies 
the critical shaking amplitude, which is analogous with 
the ID case m- There is a notable difference between 
our model in the 2D case with the model in the ID case 
[TO] . There are various shaking types. For example, a lat¬ 
tice can be shaken along one diagonal of the lattice (linear 
shaking), or elliptically (elliptical shaking), or circularly 
(circular shaking). Since separability of the system ex¬ 
ists along two primitive vectors, quasienergy dispersion 
is independent of shaking types. However, linear shak¬ 
ing preserves time-reversal (TR) symmetry, while ellipti¬ 
cal or circular shaking breaks TR symmetry. Analogous 
to orbital Hund’s rule m , there is the largest interac¬ 
tion energy at fixed momentum for repulsive bosons with 
linear shaking than for other shaking types. Together 
with inhomogeneous band mixing, we predict the small¬ 
est critical shaking amplitude for linear shaking. Then we 
construct a low-energy effective theory to describe phase 
transitions. A critical correlation length exponent is cal¬ 
culated by the momentum shell renormalization-group 
(RG) method. In the end, we study the existence of a 
Bose-Einstein condensate (BEC) in a general shaken lat¬ 
tice system. A lot of effort has been devoted to realizing 
quantum states which are not Bose condensed mm- 


2 


We find a Bose liquid with an anisotropically algebraic 
order in a three-dimensional lattice with two directions 
shaken, and we propose to change the BEC into a non- 
condensed Bose liquid via tuning shaking amplitude ap¬ 
proaching the critical value. 

The paper is organized as follows. In Sec. [TTJ we in¬ 
troduce the model for bosons in a shaken optical lattice. 
In Sec. |III| we calculate the quasienergy spectrum and 
obtain a finite-momentum superfluid phase. Next, we 
study the interaction effect on this phase in Sec. IV and 
|Vl In Sec. |VI[ we construct a low-energy effective field 
theory to study the quantum criticality of the phase tran¬ 
sition. The existence of a BEC in a general shaken lattice 
system is discussed in Sec. m Finally, conclusions are 


presented in Sec. VIII 



II. MODEL 

The system we consider is two counter-propagating 
laser beams along the x direction and two along the y 
direction, which forms a square lattice. The lattice is 
shaken by time-periodically modulating relative phase 
<p x (t) between laser beams along the x direction and 
<p y (t) between that along the y direction via acousto¬ 
optic modulators. The Hamiltonian reads 

H{t) = ^+Vcos 2 ik r x+^-) + Vcos 2 ik r y+^^), 
z m z z 

(1) 

where k r is photon momentum, ip x (t) = f cosc^t, <p y (t) = 
f cos(u;t -by?), / is the shaking amplitude, and p is the 
relative phase between ip x and p y . <p = 0 or tt means pre¬ 
serving TR symmetry, while <p ^ 0 and 7 r means break¬ 
ing TR symmetry. A = f / (2k r ) is the maximum lattice 
displacement along the x or y direction. This model is 
separable along the x and y direction. 

Taking a transformation x x — Acoscjt, y y — 
A cos (pot + <£>), the Hamiltonian in the comoving frame 
reads 

H(t,p) = -^—p 2 + V cos 2 (k r x) + V cos 2 (k r y) — ——, 

2m m 

( 2 ) 


where the effective vector potential is A(t) = 
rmoA(sincot, sm(cot + p)). Neutral particles will act as 
charged particles in a static square lattice and an ac elec¬ 
tric field E = — moo 2 A (cosoot, cos(oot + p)). The effective 
charge is set to be unity. 

The first three static terms in Eq.([2) give a static band 
structure c\(k) and corresponding Bloch wave function 
which will serve as basis in the following anal¬ 
ysis. In this paper, we consider shaking frequency uj is 
a little blue-detuned from p x and p y bands. Moreover, 
we notice higher bands couples with 5-band via higher- 
order processes due to symmetry and hence only keep 5, 
p x and p y bands. And we numerically verify that our fol¬ 
lowing qualitative results do not change when counting 


FIG. 1: Shaking induced couplings in real space. The circles 
(yellow) denote s-orbital, the colored shapes laid along the 
horizontal direction and that laid along the vertical direction 
denote the p x and p y orbitals, respectively, and + and — signs 
denote the orbital phase. 


higher bands. In these bases, the tight-binding form of 
the Hamiltonian in the comoving frame is given by 

fc 

hA 

F,fc 

( 3 ) 



where ' 1 ^ and are creation and annihilation op¬ 

erators of a particle with quasimomentum k in A band, 
respectively, and A is p x ,p y or s. 

The Hamiltonian in momentum space is given by 


Hk(t, <p) = 



where 

tlsp T COS k x ,yj 

h s = (w StPy (r)\ip x \w S ' Py (r - e x )) 

h P = ( w pA r )\iPx\w Px (r ~ e x )), 
h sp = {w Px (r)\ip x \w s (r)), 
h sp i = {w Px (r)\ip x \w s (r - e x )), 


x 


( 4 ) 


( 5 ) 

(6) 

( 7 ) 

(8) 
( 9 ) 


where Ay(f) is the 7 component of the effective vector 
potential A(t), e 7 is the primitive vector along the 7 di¬ 
rection, 7 is x or y, w\(r) is the Wannier function of the 
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FIG. 2: Quasienergy spectrum with V = 13 E r ,huj m 

6.6E r ,p = 0, where E r — Ti 2 k 2 /(2m) is the lattice recoil 
energy, (a) Quasienergy spectrum before shaking. The half 
transparent surface denotes the dressed s band with energy 
lifted by Tiuj. (b,c) Quasienergy dispersion of the uppermost 
band for (b) / = 0.02 and (c) / = 0.08. 


A band, A is p x , p y or 8, and (• • • | • • • | • • •) denotes inte¬ 
gral in the coordinate space f dr • • •. Real coupling am¬ 
plitudes h s ,h p , h sp i, h sp denote shaking induced nearest- 
neighbor hopping between s bands, between p bands, be¬ 
tween s and p bands, and onsite coupling between s and p 
bands, respectively, as shown in Fig. [l] The first matrix 
in Eq.Q represents a static band structure. And the last 
two matrices represent shaking-induced coupling among 
three bands. It is essential that lattice shaking induces 
hopping between s and p bands, which is symmetry for¬ 
bidden in the absence of shaking. Here shaking plays the 
role of external field breaking inversion symmetry, which 
is similar to mixing the p z band with p x , y band by an 
external electric field in the orbital Rashba effect [28 . 


III. FINITE-MOMENTUM PHASE 

By diagonalizing the Floquet operator, i.e., time- 
revolution operator in a time period Tq = 27t/cj, 

U(T 0 , <p) = te^-fo 0 ( 10 ) 

one can obtain a quasienergy spectrum as shown in Fig. 
[ 2 ] There is a certain critical shaking amplitude such 
that for /</[?, the uppermost band exhibits a single 
minimum at zero momentum, and for />/]?, the upper¬ 
most band exhibits four minima at finite momenta. 

Symmetry of a periodically driven system must be con¬ 
sidered at the Floquet operator [/(To, p) level [29 . Since 
U (Tq, p) is separable along x and y direction, quasiener¬ 
gies do not depend on the relative phase p. For p = 0, 
the original Hamiltonian H(t,p = 0) in Eq. ([ 2 ]) has D 4 
symmetry, so does the quasienergy spectrum. So the up¬ 
permost band dispersion has D 4 symmetry for any p. 

To describe the system, an effectively static Hamilto¬ 


nian H e ff is defined as 

u(T 0 ,(fi) = (n) 

We will analyze a rotating-wave-approximation (RWA) 
Hamiltonian, i.e., the zero order (1 /oj)° term of the ef¬ 
fective Hamiltonian [30] . which is given by 

H£ WA [}P) = (0\t)(H k (t)-id t )0(t)f } 

e Px (k) 0 Cl(k x ) \ 

_ 0 e P Jfc) e-^Q(ky) , (12) 

£l(k x ) e l(p Q(k y ) e s (k) + uj ) 

where 


&{kx,y) — 2 ^i^x,y) 5 

(13) 

/ 1 0 0 \ 


0 (t) =01 0 

(14) 

\00 e iult ) 



and superscript (0) denotes the static part. Coupling 
strength £l(k x , y ) is proportional to shaking frequency and 
amplitude. The RWA Hamiltonian in Eq. (p~2]) indicates 
that lattice-shaking-induced couplings result in level re¬ 
pulsion effect, which is the strongest along k x = 0 and 
k y = 0 directions. The level repulsion effect combined 
with T >4 symmetry will give rise to four global minima 
at finite momenta (±fc c , ±fc c ) instead of one at zero mo¬ 
mentum in the uppermost band as shaking amplitude / 
increases. 

The finite-momentum BEC has been proposed in a 
non-separable square lattice subjected to the off-resonant 
shaking [31]. Our model has more orbital physics, which 
will be shown in SecIVl 

IV. SPONTANEOUS SYMMETRY BREAKING 

Let us consider interacting bosons condensing at the 
finite-momentum state with minimal kinetic energy in 
the uppermost band. The interaction reads 

Hintit) = g I dr&{r,t)&(r,t)9(r,t)9{r,t), (15) 

where \ht(r,£) and T(r,£) are creation and annihilation 
operators of the condensate state, respectively, and pos¬ 
itive g is the repulsive interaction strength. 

Bosons can either condense at one of the four degener¬ 
ate finite-momentum states or at the superposition state. 
Assume the single-particle ground state is a superposition 
state 

4 

V’fecC 7 *^)( 16 ) 

i= 1 

where r ipk i (r,t) is one of the four degenerate states, ki = 
(±fc c , ±fc c ) is the condensate momentum, and constant 
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cti satisfies J2i=i \ a i\ 2 = 1- One can write i/j ki (r,t) in the 
comoving frame as 

Tpki(r,i) = b Pxikc 4> Pxik .(r) + b Py , kc (l> Pyiki (r) 

+ e LUJt b s , kc (j) s ,k i (r), (17) 

where b\, kc is a combination coefficient and dependent 
on A and fc c , and A is p x ,P y , or s. 

The time-average mean-filed interaction energy per 
particle condensing at the superposition state i/j kc (r,t) 
in the laboratory frame is given by 

tint(k c ) = ^J^ dtg J dr\-ip kc {r,t)\ 4 . (18) 

By minimizing the interaction energy ei nt (k c ) with 
respect to {a*}| =1 , one obtains ( 01 , 02 , 03 , 04 ) = 

(±1,0, 0,0) or (0, ±1,0,0) or (0,0, ±1,0) or (0,0,0, ±1). 
Since there is only one nonzero cq, we neglect the phase 
of a,i. So bosons only condense at one of the four finite- 
momentum states, which breaks D 4 symmetry sponta¬ 
neously. 

In the process of turning on shaking adiabatically, 
bosons will remain in the uppermost band. When shak¬ 
ing amplitude across a critical value, phase transition 
from the NSF phase to the D 4 SF phase happens. 


V. CRITICAL SHAKING AMPLITUDE 

The interaction effect also modifies the critical shaking 
amplitude. By minimizing the total energy that consists 
of kinetic and interaction energies of bosons in the up¬ 
permost band with respect to quasimomentum, one can 
obtain the condensate momentum (±fc c , ±fc c ). Here the 
methods we use to calculate kinetic and interaction en¬ 
ergies are the same as the methods used in Sec. |III| and 
[IV| respectively. When k c turns out to be nonvanishing 
with the increasing shaking amplitude, critical shaking 
amplitude f c is obtained. 

Fig. |3ja,b) show repulsive interaction effect enlarges 
the critical shaking amplitude f c in the deep lattice limit. 
This is because inhomogeneous band mixing in momen¬ 
tum space causes a global minimum of interaction energy 
at zero momentum in the deep lattice as shown in Fig. 
[3|c), which has been illustrated in the ID case [TO]. In¬ 
stead, a local maximum has also been predicted at zero 
momentum in the shallow lattice limit, which leads to a 
smaller critical shaking amplitude f c < [10]. We will 

focus on the deep lattice case in this paper. 

Besides, Fig. [3^a,b) also show shaking types can mod¬ 
ify critical shaking amplitude f c . The k c — f curve for 
the interaction case in Fig. [3ja) and the gn — f c curve in 
Fig. |3jb) change with relative phase p and are bounded 
by corresponding curves with p = 0 or tt and that with 
(p = ± 7 r/ 2 . 

One can write the eigenstate of the uppermost band in 
the comoving frame as 

Tpk(r,t,ip) = Cp x ,k<t>p x ,k{ r ) + e ~ l<p c Pyt k(f ) Py,k(r) 



FIG. 3: (a,b) Interaction shifts of critical shaking amplitude 
with V = 13 E r ,uj = 6AE r . The solid line denotes p = 0 
or 7 r, and dashed line denotes p = ± 7 r/ 2 . (a) Condensate 

momentum component k c as a function of shaking amplitude 
/ for non-interacting (blue), and interacting (red) cases with 
gn = 0.04E r . (b) Phase diagram for a given shaking fre¬ 

quency u. Left region: NSF phase. Right region: D 4 SF 
phase, (c) Interaction energy {tint) with p = 0 and ( f,gn ) 
marked as a star in (b). 


±^ Cs^k^s^k {?), 


(19) 


where cy*, is the combination coefficient of the eigenstate 
(p = 0) in the uppermost band and A is p x ,p y 


of H* wa 


or s. Since Hj? WA ((p = 0) is real and symmetric, c\, k 
can be set to be real. Here we use the RWA Hamiltonian 
for simplicity of analysis. 

Using Eq. (18), the difference between time-average 
interaction energy per particle in a system with broken 
TR symmetry and that in a system with TR symmetry 
is given by 


Ae^ n t(^5 P) £int(ki P — 0) 

= -4^sinV4 )fe y )fe [/^, (20) 

where v is the site occupation number and JJ^ xPy = 

g J drU^, k {r)(pp y ,k{ r )\ 2 - 

Eq.(20) shows two important ingredients. One is 
— sin 2 cp signifying level of TR symmetry breaking. 

Generally, when one Hamiltonian with TR symmetry is 
unitarily transformed to another Hamiltonian, the spec¬ 
trum rather than TR symmetry is always invariant. As¬ 
suming H = HQ with TR symmetric Hamiltonian H 
and unitary operator Q, one obtains 

[T, H] = (Q f HQ - Q ] HQ)T, (21) 


where T is the TR operator and Q = TQT x . If Q 

commutes with T, then Q = Q and H has TR symmetry. 
Otherwise TR symmetry is broken generally. In our case, 
H ~ Hj} WA ( 0) has TR symmetry, Q has the matrix form 


/ 1 0 0 
0 0 
\ 0 0 1 


( 22 ) 
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Quasi-momentum k /k 

x r 

FIG. 4: Contour plot of Cp xk Cp yk U k x,Py , the momentum- 
dependent part of the interaction energy difference 
A ei n t(k, p), with parameters same as used in Fig. [3^c) . 


and commutation relation [T, H k WA (ip)\ is proportional 
to sin ip. So TR symmetry is conserved for ip = 0 or 
7r and broken maximumly for ip = ±7t/2. A €i nt (k,ip) 
in Eq. (20) at the fixed momentum decrease as the level 
of TR symmetry breaking increases, which is similar to 
orbital Hund’s rule [15]. 

The other key ingredient of interaction energy differ- 


Px iVy 


IS 


ence in Eq. (20) is momentum dependence. U t 
a positive constant in the deep lattice limit. Inhomo¬ 
geneous band mixing in momentum space causes the 
fact that the momentum-dependent part Cp x ^Cp y ^U k x ' Py 
has a global maximum at zero momentum, as shown in 
Fig.@. Here we can see ip appears only in interaction 
components involving p x and p y orbitals. The reason 
is that other terms involving cp, such as c^ k c 2 yk e ±l2ip , 


contain factor e ±l2ujt because of the energy difference be¬ 
tween s and p bands and can be neglected in the sense 
of time average. 

The two ingredients determine curvature of interaction 
energy at zero momentum is minimum when cp = 0 or tt 
and maximum when ip = ± 7 t/ 2 . So there is the smallest 
critical shaking amplitude for ip = 0 or tt and the largest 
one for ip = ±7t/2. 


VI. EFFECTIVE FILED THEORY 

In this section, we introduce a low-energy effective ac¬ 
tion to describe all three phases, i.e., NSF, D 4 SF and MI 
phases. Based on this action, we will show phase diagram 
calculated by mean-field theory and a critical exponent 
calculated by momentum shell RG theory. 

In order to construct the effective action, two impor¬ 
tant factors from microscopic analysis above must be 
considered. First, kinetic energy has the quartic form 
of k\ + ky + a(k 2 + k 2 ) at small momentum. Second, 


momentum-dependent interaction has the quadratic form 
of a + f3(k 2 + k 2 ) at small momentum. These two factors 
also agree with D 4 symmetry. The low-energy effective 
action of the d-dimensional lattice shaken along the x 
and y directions can be written as 

nl/ T n 

£($,$*)= / dr / & d r{K l ^*d T ^+K 2 \d T ^\ 2 +£{^^*)}, 

(23) 

where 


£($,$*) = |d 2 $| 2 + |d 2 $| 2 + a|V $| 2 + T + r |$| 2 

+ a|$| 4 +/3|$V$| 2 , (24) 


T is temperature, is superfluid order parameter, V = 
(d x ,d y ), T = 0 for d = 2 and T = \d z $\ 2 for d = 3. The 
signs of parameters a and r can be inverted by tuning 
shaking amplitude / and interaction strength g , respec¬ 
tively. Parameter a is considered to be positive for repul¬ 
sive interactions. And parameter /3 can be either positive 
in the deep lattice limit, or negative in the shallow lattice 
limit. We will take /3 > 0 for simplicity. 

Assuming = |ff>|e zfc r , £ can be rewritten as 

£(|$| ,k) = (fc 4 + fc 4 + afc 2 +r)|$| 2 + (a + /3fc 2 )|$| 4 , (25) 


where k = ( k x , k y , • • •) and r = (x, y, • • •) are d- 
dimensional vectors. By minimizing £ with respect to 
|<F| and fc, one obtains three different phases: (1) the MI 
phase with = 0; (2) the NSF phase with ^ 0 and 
k = 0; (3) the D 4 SF phase with <F / 0, k x / 0, k y / 0 
and k 2 — k 2 —k 2 = 0. Phase boundaries are also obtained: 

(1) r = 0 and a > 0 separating the NSF and MI phase; 

( 2 ) r = 2 aa//3 < 0 separating the NSF and D 4 SF phase; 

(3) r = a 2 / 2 and a < 0 separating the D 4 SF and MI 
phase, which is different from the ID case m- There is 
a mean-field tricritical point (a,r) = (0,0). In the vicin¬ 
ity of this tricritial point, a and r are proportional to 
f — fc and g~g c , respectively. The phase diagram in /- 
and g-terms is shown in Fig. ©• 

Next we will study the critical correlation length ex¬ 
ponent within momentum shell RG approach. 

At zero temperature and in (d + 1 ) dimensional mo¬ 
mentum and frequency space, the action in Eq.(23) can 
be rewritten as 


£(*.**) = / i0p^*{^“){-iKiu + K 2 u 2 

+ T k y + u(k 2 + k 2 ) + 7/c T t to) 

+ / {<a + P(k2 X kAx + kff) 

J kco 

Q*(k2,U2Mk3,U S )$(k4,(jJ4), (26) 


where Tk = 0 for d = 2, Tk = k 2 for d = 3, f koj = 

f k 2 — k 3 — k 4 )S(uJi -h UJ 2 — Cd 3 — UJ 4 ) 

and A denotes high momentum cut-off. 
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Interaction (g-g c )/g c 

FIG. 5: Mean-field phase diagram. g c is the critical inter¬ 
action strength for the NSF-MI transition. /? is the criti¬ 
cal shaking amplitude calculated by minimizing the single¬ 
particle quasienergy dispersion. 



(a) (b) (c) 


FIG. 6: The one-loop Feynman graphs. Graph (a) con¬ 
tributes to renormalizing parameters a and r. Graph (b) and 
(c) contribute to renormalizing parameters a and /3. 


The one-loop Feynman graphs for renormalizing the 
parameters in Eq.(26) are shown in Fig. ©• The study 
of the critical exponent is divided into two cases. 

Case A. Without particle-hole symmetry. K\ ^ 0, 
so the i^-term becomes irrelevant. In this case, scaling 
dimensions of the parameters read 


[k x \ = [k y \ = I, [k z ] = 1 , [a] = 1 , [r] = 2 , [w] = 2 , 
[a]=3-d,[/3]=2-d,[$] = -^. (27) 


The upper critical dimension is 3, which is 1/2 larger than 
that in the ID shaken lattice m due to the extra shaking 
direction. For d = 2, the <a-term is relevant and the /3- 
term is marginal. It’s different from the irrelevant /3-term 
in the ID case m- So we need to consider corrections 
from the /3-term. For graphs in Figj 6 ] (b,c), we need to 
expand them in powers of external momenta. The one- 
loop RG flow equations read 


da 

d7 = a ’ 

dr 

d[ = 2r ’ 
da 

— = ea - 

1 + r 


a 


h(a), 

a/3 


h(p) + 


—a/3L 2 (a)\ - 


a 


(1 + r)* 


(1 + r ) 2 
M 2 (a), 


\a 2 J 2 (a) 


(28) 

(29) 

(30) 

) 

(31) 


where 


Ha) = [ 

Jo 

J 2 (a) = [ 

Jo 

£2 (a) = f 

Jo 

M 2 (a) = [ 
Jo 


2ir d p ~ a + \/ a2 + 4(cos 4 p + sin 4 p) 
(27r) 2 4(cos 4 p + sin 4 p) 

dp -a + yja 2 + 4(cos 4 p + sin 4 p) 


(27T) 


2 (cos 4 p + sin 4 p) 


a + 3 cos 2 p- 


a + y a 2 + 4(cos 4 (p + sin 4 p) 
cos 4 p + sin 4 p 


dp -a + Ja 2 + 4(cos 4 p + sin 4 p) 


(27T) 


2 (cos 4 p + sin 4 0 ) 


d</> -a + y a 2 + 4(cos 4 p + sin 4 0) 


(27T) 


cos 4 p + sin 


a + cos 2 p- 


a + cos 2 p- 


fl + ya 2 + 4(cos 4 p + sin 4 0) 
cos 4 p + sin 4 p 


(32) 

(33) 

cos 2 p, (34) 


a + y a 2 + 4(cos 4 p + sin 4 0) 
cos 4 p + sin 4 p 


cos" 


(35) 


and e = 3 — d. 


flow equations are given by 


The nontrivial fixed point lies at (r*, a*, a*, /3*) = 
(0,0, 7 ^ 0 ) J 2 ^/^/q) 2 2 ^ g2 )• Defining new variables = 
r—r*,5a = a—a* ,5a = a—a*, 5/3 m /3—/3*, the linearized 


/ £r 
d £a 
dZ ^ 

W/3 


/ 2 0 0 

0 1 0 

0 0 -e 


\0 0 2 


j 2 (o)-m 2 (o)^ 
J 2 (o) e 


-1 - 


J 2 (o) e / 


/ \ 

5a 

5a 

\6/3J 
(36) 
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Eigenvalues of the matrix in Eq.(36) are 2 , 1 , —e, —1 — The upper critical dimension is 4. For d — 3, the one-loop 
eZ/ 2 ( 0 )// 2 ( 0 ). Then the scaling dimension of parameter RG equations read 
r at the nontrivial fixed point is y r = 2. The correlation 
length exponent of the superfluid transition can be calcu¬ 
lated as v = l/y r = 1 / 2 . It is the same as the mean-field 
value in the usual Bose gas |32, 33. and the value in the 
ID case with K\ ^ o m because of no one-loop correc- 


da 


p 


— = a -f ,- 

dl 2y/l -Rr 


h{a ), 


tions on r from interaction as shown in the flow diagram _ — -\ _ - 


m Figj7^a). dl 

Case B. With particle-hole symmetry. K\ = 0. In this q a 

case, scaling dimensions of the parameters read 


[kx\ = [ky\ = X, [k z \ = 1, [a] = 1, [r] = 2 , [u] = 1, 


M = 4 - d, [/3] = 3 - d, [4>] = - 


d T 2 


(37) 


vT+r 

1 

— = ea -; 

d 2 (1 + r)i 


(3 


2alz(a) + —Jz^a) 


(38) 

(39) 


-a 2 I^{a) + aPJz(a) + fi 2 Lz(a) 


(40) 


^ = (e - 1)0 - - - T 

dl 2(1 + r) 2 


OLfih{a) + yJ 3 (a) 


+ 


(1 + r)a { 2 


-<a 2 [a/ 3 (a) + 3J 3 (a)] + a/3M 3 (a) + 0 2 7V 3 (a) 


+ 


1 


(1 + r) 


y [<a 2 /M a ) + Qif3R 3 (a) + / d 2 5' 3 (a)] , 


(41) 


where 


sind, 


I sin 2 9 


+ cos 2 d 


J3(a) = (S? l M l yy+cos 2 " 


sin 2 d 


— (a sin 2 d + cos 2 d) + y (a sin 2 d + cos 2 d ) 2 + 4 sin 4 d(sin 4 0 + cos 4 0) 

2 sin 4 d(sin 4 0 + cos 4 0 ) 

(42) 

— (a sin 2 d + cos 2 d) + (a sin 2 d + cos 2 d ) 2 + 4 sin 4 d(sin 4 0 + cos 4 0) 


2 sin 4 d(sin 4 0 + cos 4 0 ) 


sin 3 d, 


(43) 


-| p7T pZ 7T 

sin 5 d, 


I sin 2 d 


+ cos 2 d 


— (a sin 2 d + cos 2 d) + y (a sin 2 d + cos 2 d ) 2 + 4 sin 4 d(sin 4 0 + cos 4 0) 


2 sin 4 d(sin 4 0 + cos 4 0 ) 


(44) 


Ms (a) = 


pTT p2lT 

l M l 


I sin 2 d 


+ cos 2 d 


a(l + cos 2 0 ) + cos 2 0(7 + cos 20 ) 


sin 3 d 


— (a sin 2 d + cos 2 d) + y (a sin 2 d + cos 2 d ) 2 + 4 sin 4 d(sin 4 0 + cos 4 0) 

2 sin 4 d(sin 4 0 + cos 4 0 ) 

— (a sin 2 d + cos 2 d) + J (a sin 2 d + cos 2 d ) 2 + 4 sin 4 d(sin 4 0 + cos 4 0) 


(45) 


2 sin 2 d(sin 4 0 + cos 4 0 ) 


^Vs(a) = 


16(2t r) : 


/>7T p2lT 

l de l d *\ 


I sin 2 d 


+ cos 2 d 


sin 5 d 


a(3 + 2 cos 20) + cos 2 0(5 + 2 cos 20) 


— (a sin 2 d + cos 2 d) + y (a sin 2 d + cos 2 d ) 2 + 4 sin 4 d(sin 4 0 + cos 4 0) 

2 sin 4 d(sin 4 0 + cos 4 0 ) 

— (a sin 2 d + cos 2 d) + a/( a sin 2 d + cos 2 d ) 2 + 4 sin 4 d(sin 4 0 + cos 4 0) 


sin 2 d(sin 4 0 + cos 4 0 ) 


( 46 ) 






























































PM = 




2?r L m 2 0 


(27r) 3 
sin 3 9 cos 2 0 


d<M 


s 2 0 


(a sin 2 6 + cos 2 0) + \ (a sin 2 0 + cos 2 0) 2 + 4 sin 4 0(sin 4 0 + cos 4 0) 


a + cos 


i? 3 (a) = 


(27t) 3 
sin 5 0 cos 2 0 


/*7T /»27T 

l de l ^ 


I sin 2 0 


+ cos 2 0 


a + cos 


5 3 (a) = 


8(2tt) 3 
sin 7 0 cos 2 (j) 


P7 T /‘2tT 

l M l d *1 


I sin 2 0 


cos 2 0 


2 sin 4 0(sin 4 + cos 4 (j>) 


2 ~\~ \/( a s ^ n2 ^ + cos2 ^) 2 + ^ s ^ n4 $( s i n4 0 + cos 4 <j)) 

> (47) 

sin 2 0(sin 4 </> + cos 4 <fi) 

— (a sin 2 0 + cos 2 0) + yj (a sin 2 0 + cos 2 0) 2 + 4 sin 4 0(sin 4 (j) + cos 4 0) 

2 sin 4 0(sin 4 0 + cos 4 <fi) 

2 

2 6) + ^/(a sin 2 0 + cos 2 6) 2 + 4 sin 4 0(sin 4 </> + cos 4 0) 

, (48) 

sin 2 0(sin 4 </> + cos 4 <fi) 

— (a sin 2 0 + cos 2 9) + yj (a sin 2 0 + cos 2 0) 2 + 4 sin 4 0(sin 4 </> + cos 4 0) 


2 sin 4 0(sin 4 (j) + cos 4 0) 


a + cos 2 <fi- 


— (a sin 2 0 + cos 2 0) + y (a sin 2 0 + cos 2 0) 2 + 4 sin 4 0(sin 4 </> + cos 4 0) 
sin 2 0(sin 4 </> + cos 4 <fi) 


(49) 


and e = 4 — d = 1. The nontrivial fixed point lies at 


(r*,a*,at*,P*) 


( 2e 27J 3 (0) + 2P 3 (0) 2 2e 
V 5’ 257 3 (0) ’ 57 3 (0) ’ 

2[27J 3 (0) + 2P 3 (0)] 

25/1(0) ) ■ (5C 


Defining Sr = r — r*,£a = a — a*,£a; = a — a* and 
S/3 = (3 — f3* gives the linearized equations 


/ Sr 
d Sa 
d l Sa 

\sp 


/2-§e 
0 
0 

0 


5/3 (0) 


3 / 3 ( 0 ) 1 


<9a 

1 

0 


i=oe 


2/ 3 (0)(l + |e) 
0 


—e 


[27J 3 (0) + 2P 3 (0)] 


5/3(0) c 


nPa + W 

M (1 + 


-1 


2Js(0) 

5/3(0) 6 

, M 3 (0)+/? 3 (0) 
+ / 3 ( 0 ) 



(51) 


The eigenvalues of the matrix in Eq.([5l)) are 2 — 
2e/5,1, -e and -1 + 4e/5 + 2e(M 3 (0) + P^0))/5/ 3 (0). 
Here it is the important difference from Case A that the 
r-term gets corrections from interaction as shown in flow 
diagram in Fig. ([7]) (b). The correlation length expo¬ 
nent of the superfluid transition is v = 1/(2 — 2e/5) = 
1/2 H- e/10 = 3/5. It is different from bosons with quar- 
tic dispersion in only one direction with particle-hole 
symmetry m due to different upper critical dimensions, 
while it is the same as conventional bosons with quadratic 
dispersion with particle-hole symmetry in three dimen¬ 
sion, and thus, belongs to the 0(2) rotor model class, up 
to e order [34] . 

Systems with parameters r and a that lie in the right 
region of the critical surface in Fig. 0 will eventually 




FIG. 7: RG flow diagrams of (a) case A and (b) case B on 
the r — a plane. 


flow towards the MI phase, whereas systems with r and a 
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in the left region of the critical surface will flow towards 
the NSF phase for the initial a > a* or the D 4 SF phase 
for the initial a < a*. 


VII. THE FATE OF THE BEC IN GENERAL 
SHAKEN LATTICE SYSTEMS 

In sections above, we suppose a BEC exists in the sys¬ 
tem. However, it is known that low-energy density of 
states can be increased by SO coupling [271I35H38] or lat¬ 
tice shaking [27, . Fluctuations may destroy off-diagonal 
long-range order (ODLRO). In this section we will first 
study existence of ODLRO in the 2 D cases, and then in 
general cases. 

Let us first assume bosons condensing in the NSF 
phase and write <f> = po + Spe 10 , where po = —r/2a 
is superfluid density, Sp is density fluctuation and 0 is 
phase fluctuation. Then we substitute the <f> field into 


Eq. (23) and expand the action in Eq. (23) to quadratic 
order in Sp and 0. By integrating out Sp field, the low- 
energy effective action for 6 field is given by 

S eff (9) =p 0 j d d rdr{K(d T e) 2 + d$6 + d d 0 + a(W) 2 
+ t|, (52) 

~ j ^-2 

where K = ^ + K 2 , a = a + 0p o - 

In Gaussian approximation, the correlation function 
can be written as 

($*(r)$(0)) =poe“5 <(e(r) “ 0(o))2) . (53) 


At finite temperature, Eq. (53) can be written as 


($*(r)$(0)) = po ex P< — T 


■/ 


d d k 


(1 - e lkr ) 


{2it) d 

+ ky + a(kl + ky ) + Tk) j • (54) 

For d = 2, in the large separation limit, the integral in 
Eq. (54) can be approximated by 


n2ir />/ 

/ d ° 

J 0 Jo 


dk 


1 


(27r) 2 k 3 ( cos 4 0 + sin 4 6) + dk 


(55) 


As k —» 0, the integrand in Eq. ([55]) behaves as 1 /k for 
a > 0 and 1/k 3 for a = 0 . So the integral is divergent, the 
correlation function in Eq. (54) approaches zero at large 


separation, and ODLRO is absent at any finite tempera¬ 
ture. 

For d = 3 and a > 0, in the large separation limit, the 
integral in Eq. (54) can be approximated by 


/ 


crfc 


( 2 tt ) 3 k* + k* + a(kl + fc|) + k \' 


(56) 


After integrating over k z , the integral above is given by 


/ d9 

Jo Jo 


dk 


Sir 2 


Y/ k 2 ( cos 4 0 + sin 4 0) + a 


(57) 


As k 0 , the integrand in Eq. (57) behaves like 1 . So 


the integral is finite and ODLRO exists. 

For d = 3 and a = 0, the integral in the exponent in 


Eq. (54) reads 


/ 


d 3 fc 


ikr 


(27r) 3 K 


K 


k 2 


(58) 


The critical point a = 0 means critical shaking amplitude 
f c in Sec. [V| i.e., phase boundary separating NSF and 
D 4 SF phases in the mean field level in Sec. |VI| At large 
separation, the correlation function is given by 


(4>*(x)4>(0)) - 

<<l>*(*)$(0)> ~ k|-" T . 


(59) 

(60) 


where p = r (5/4) 2 /7 r 5 / 2 and F(z) denotes Euler gamma 
function. So there exists BEC only at zero tempera¬ 
ture and noncondensed Bose liquid at finite temperature. 
This algebraically ordered Bose liquid is anisotropic. 
Since fluctuation is enhanced by shaking, the correlation 
function decays faster along the shaking directions. 

At zero temperature, Eq.( 53) can be written as 


($*(r)$( 0 )) = poexps - 


W 


d d kdoj 


y (1 — e ikr ) 


(27r) d + 


^Kcd^ + A: 4 + ky + a(k 2 + ky) + 7 


-1 


(61) 


The zero temperature results are equivalent to adding a 
unshaken direction to the corresponding finite tempera¬ 
ture case. The correlation function in Eq.(61) remains 
finite at large separation for a > 0 and d = 2 and 3. At 
the critical point a = 0 and in the large separation limit, 
the vanishing correlation function is found for d = 2 , 
which is consistent with results in SO-coupled BEC with 
similar dispersion EH HU, whereas the finite correlation 
function is found for d = 3. 

From calculations above, we know phase fluctuation 
destroys ODLRO at any finite temperature for d = 2 and 
a > 0. The effect is even stronger at the critical point due 
to pure quartic dispersion. For d = 2 and a = 0, ODLRO 
does not exists even at zero temperature pa mi. For a 
system with d = 3, there exists ODLRO when a > 0 
and quasi-long-range order at finite temperature when 
a = 0. Therefore the BEC can be changed into a non¬ 
condensed Bose liquid by tuning the shaking amplitude 
/ approaching the critical value f c . 

Raman-induced SO coupling and lattice shaking have 
generated quartic dispersion [27] . And higher-order 
terms may be generated in the future. Next we will study 
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TABLE I: Existence of ODLRO for a = 0 at finite tempera¬ 
tures. 



Tk 

K + % 
(n > 4) 

K + K + Tu 
(m > n > 4) 

k” + k™ + k l z 
(l > m > n > 4) 

d=2 

N 

N 

N 

— 

CO 

II 

"O 

O 

O 

N 

N 

TABLE II: 

ture. 

Existence of ODLRO for a - 

= 0 at zero tempera- 


T k 

K + Tk K + k™ + % 

(n > 4) ( m > n > 4) 

K + k™ + k l z 
(l >m>n> 4) 

d= 2 

0 

0 

N 

— 

d= 3 

0 

0 

O: n — m — 4 

O or n = 4, m = 6, / < 10; 

N: otherwise 


the feasibility of changing the BEC into a non-condensed 
Bose liquid in a system with a general dispersion. 

The existence of ODLRO is obtained by checking if the 
correlation function in Eq. ( [53] ) is finite in large separation 
limit. And the results at the critical point a = 0 are 
shown in Table [T] and [TTJ where N represents having no 
ODLRO and O represents having ODLRO. For a > 0 , at 
finite temperatures, ODLRO exists only in systems with 
d = 3. So systems with dispersion k^^k^ -\-a{k 2 3 4 5 6 7 8 9 10 11 12 -\-ky)-\-k 2 


or k™ + k™ + k l z + a{k 2 + k 2 + k 2 ) with l > m > n > 4 can 
be used to change the BEC into a noncondensed Bose 
liquid by tuning the shaking amplitude approaching the 
critical value f c . 


VIII. COCLUSIONS 

In conclusion, we have investigated quantum phase 
transition of bosons in a shaken lattice by using Floquet 
theory and low-energy effective field theory. We found 
there was a D 4 SF phase with spontaneous symme¬ 
try breaking and calculated the critical shaking ampli¬ 
tude f c for the NSF-D4SF phase transition. We further 
demonstrated both the interaction effect induced by in¬ 
homogeneous band mixing and the shaking types could 
modify f c . We identified a quantum tricritical point of 
NSF, D 4 SF and MI phases and studied quantum critical¬ 
ity nearby the tricritical point. And the critical exponent 
is expected to be measured by in situ density measure¬ 
ments 23 in the future. Moreover, we found anisotropi- 
cally algebraic order and proposed to turn the BEC into 
a noncondensed Bose liquid by tuning the shaking am¬ 
plitude approaching the critical value f c . 
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